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SUMMARY 


Two  valve  test  systems  were  studied.  They  each  use  a  stored  energy 
approach;  one  is  rotational,  the  other  is  translational  (drop).  An  idealized 
dynamic  analysis  indicates  the  rotational  system  will  generate  peak  forces  and 
torques  that  will  require  a  massive  test  fixture.  The  test  fixture  needs  a 
thorough  design  and  dynamic  analysis  to  assure  it  will  operate  under  such 
severe  dynamic  forces. 

The  drop  test  will  almost  equal  the  performance  of  the  rotational 
system.  It  also  appears  to  be  less  costly.  However,  it  is  near  the  practical 
limits  of  ball  weight,  drop  height  and  available  reaction  structure. 

Neither  system  produces  the  design  discharge  for  enough  time  to  evaluate 
the  regulating  valve's  performance.  Precise  timing  of  the  start  of  pumping  in 
the  rotational  test  is  critical  to  attaining  the  desired  discharges. 


PREFACE 

This  study  was  conducted  in  January  1986,  by  personnel  of  the  US  Army 
Engineer  Waterways  Experiment  Station  (WES),  under  the  sponsorship  of  the 
Defense  Nuclear  Agency  (DNA)  in  support  of  the  Silo  Test  Program-Shock  Iso¬ 
lation  Systems.  The  DNA  project  officer  was  Mr.  James  Cooper.  Mr.  Larry 
Selzer,  Aerospace  Corporation,  proposed  the  concept  as  a  means  of  evaluatiing 
a  valve  for  a  full-scale  shock  isolation  system. 

The  investigation  was  conducted  under  the  supervision  of  Messrs.  Bryant 
Mather,  Chief,  Structures  Laboratory  (SL);  James  T.  Ballard,  Assistant  Chief, 
SL;  Dr.  Jimmy  P.  Balsara,  Chief,  Structural  Mechanics  Division  (SMD),  SL;  and 
Mr.  Robert  E.  Walker,  Project  Manager,  SMD.  This  report  was  prepared  by  Mr. 
James  B.  Cheek,  SMD. 

The  Director  of  WES  during  the  investigation  and  preparation  of  this 
report  was  COL  Allen  F.  Grum.  The  Technical  Director  was  Dr.  Robert  W. 
Whalin. 
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STUDY  OF  STORED  ENERGY  SYSTEMS 


PROPOSED  FOR  TESTING 
A  PRESSURE  REGULATING  VALVE 


PART  I:  INTRODUCTION 

1 .  The  analysis  documented  in  this  report  was  accomplished  as  part  of  a 
feasibility  study  to  develop  design  requirements  for  a  device  proposed  for 
testing  a  hydraulic  pressure  regulating  valve.  Because  cf  the  valve's  high 
operating  pressure  and  discharge,  designing  the  test  device  presents  many 
difficult  analysis  problems.  This  study  looks  at  but  one  of  those  problems  in 
a  highly  idealized  fashion.  Nevertheless,  the  analysis  is  useful  in  that  it 
establishes  the  best  performance  attainable  from  a  "perfect"  system.  That 
performance  can  be  used  to  see  how  well  it  meets  test  system  requirements. 

From  that  evaluation,  decisions  on  any  changes  needed  to  a  practical  system 
can  be  made  and  more  extensive  engineering  analysis  can  be  conducted. 

PART  II:  IDEALIZED  ANALYSIS,  ROTATIONAL 

Valve  Specification 

2.  The  following  analysis  was  done  on  a  fixture  proposed  to  test  a 
regulating  valve  at  a  constant  regulating  pressure  of  ^,350  pounds  per  square 
inch  (psi),  at  a  design  maximum  flow  of  26.^  gallons  per  second  (gps). 

Test  Fixture  Data 

3.  The  Test  Fixture  consists  of  two,  30  inch  diameter,  solid  disc 
flywheels  weighting  1,000  pounds  (lb).  Each  is  connected  to  a  crankshaft 
having  a  one-inch  offset  (throw)  which  is  in  turn  linked  by  a  connecting  rod 
to  a  pump  as  shown  in  Figure  1  below.  The  pump  cylinder's  ID  is  16  inches  and 
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Dynamic  Analysis 

4.  The  analysis  is  based  on  a  perfect  system  (no  mechanical  or  pumping 
loss).  The  results,  presented  in  Appendix  A,  indicate  the  following: 

_a.  At  424  revolutions  per  minute  (rpm),  the  pump  will  discharge 

26.4  gps,  peak ,  provided  the  pumping  starts  when  the  crankshaft 
angle  is  90°.  The  discharge  will  decline,  as  shown  in  Figure 
A2,  to  zero  in  71  milliseconds  (msec). 

b_.  Peak  torque  in  the  drive  is  near  50,000  foot-pounds  (ft-lbs) 

(see  Figures  A1  and  A3). 

c_.  The  drive  will  slow  from  424  rpm  to  zero  rpm  in  70°  of  crank 
rotation  which  is  70.8°/360°  =  .197  revolutions. 

cl.  The  force  required  to  operate  the  pump  shaft  is  almost  600,000 
lb  at  design  pressure. 

_e.  Operating  the  system  at  577  rpm  and  starting  the  test  at  zero 

degrees  crankshaft  angle  produces  the  desired  peak  discharge  and 
increases  the  time  of  regulated  discharge  (See  Figure  B7). 

5.  Appendix  B,  like  Appendix  A,  graphs  the  test  fixture  performance. 
However,  pumping  starts  at  a  crank  angle  of  zero  degrees.  Peak  pump  output  is 
only  16  gps  using  424  rpm  as  above  (see  Figures  B1  and  B2). 

6.  A  second  series  of  calculations  was  made  keeping  all  conditions  the 
same  except  the  shaft  speed  which  was  changed  to  577  rpm  in  order  to  raise  the 
peak  pump  discharge  to  near  the  design  specification  (26.4  gps).  The  results 
of  those  calculations  are  presented  on  Figures  B7  through  B12. 

7.  Those  calculations  indicate  a  longer  discharge  time  at  a  higher 
discharge.  However,  the  total  discharge  time  of  a  single  cycle  system  such  as 
t'  ,  is  unlikely  to  be  long  enough  to  thoroughly  exercise  the  test  valve. 

8.  The  computer  program  used  for  the  rotational  system  analysis  is 
presented  in  Appendix  C. 

9.  Appendix  D  presents  the  calculations  upon  which  the  dynamic 
analysis  is  based.  It  also  outlines  the  program  logic  used  in  modeling  the 
slowdown  of  the  drive  system. 
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PART  III.  IDEALIZED  ANALYSIS,  DROP  TEST 

Drop  Test  System 

10.  For  comparison  purposes,  a  swinging  ball  test  fixture  was 
evaluated.  The  system  consisted  of  a  ball  suspended  by  a  sling  from  a  fixed 
point.  When  directly  below  the  suspension  point,  the  ball  will  impact  the 
piston  of  a  hydraulic  cylinder  connected  to  the  regulating  valve.  Energy  is 
stored  in  the  ball  by  raising  it  to  height  H  has  shown  in  Figure  2. 


FIGURE  2 


Dynamic  Analysis,  Drop 

11.  The  dynamic  analysis  calculations  and  the  computer  program,  shown 
in  Appendix  E,  are  based  on  total  transfer  of  momentum  from  the  ball  to  the 
pump  system,  ie.,  the  ball  does  not  rebound.  The  results  show  the  following: 

a_.  The  drop  height  controls  the  peak  flow  rate.  Consequently,  for 
a  given  cylinder  diameter,  H  is  fixed  in  order  to  attain  the 
valve's  design  flow. 

_b.  For  a  given  cylinder  diameter  (thus  drop  height)  increasing  the 
weight  of  the  ball  increases  the  discharge  time  of  the  cylinder. 

£.  A  3.5  inch  diameter  cylinder,  a  *13.1  foot  drop  height,  and  a 
2,000  lb  ball  will  produce  the  design  peak  fl..w  followed  by  a 
linear  decline  to  zero  flow  for  78  msec. 

_d.  Figure  3  is  a  composite  plot  showing  the  discharge  versus  time 
graphs  of  the  best  rotational  test  (from  Figure  B8)  and  the  best 
drop  test.  This  plot  clearly  shows  the  superior  performance  of 
the  rotational  system. 

£.  Time-discharge  graphs  of  other  cylinder  diameters  and  ball 
weights  are  shown  in  Figure  E2. 
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FIGURE  3 


PART  IV:  CONCLUSIONS 


12.  The  following  conclusions  are  based  on  the  idealized  analysis  of 
the  two  systems. 

_a.  Provided  the  initial  rpm  of  the  rotational  test  is  increased 
from  424  to  577,  the  Drop  Test  will  do  an  inferior  job  of 
exercising  the  regulating  valve  because  of  the  Drop  Test’s 
linear  decline  in  discharge.  However,  the  drop  test  will  cost 
less . 

1).  The  drop  test  is  near  the  practical  limits  of  ball  weight  and 
drop  height.  Going  beyond  the  design  discharge  requires  a  four 
fold  increase  in  drop  height  to  double  the  peak  discharge. 
Changing  to  a  larger  cylinder  diameter  increases  discharge 
directly  with  the  square  of  D,  but  decreases  flow  time  with  D^. 

c_.  The  very  short  test  time  of  the  rotating  system  produces 

extremely  large  forces  in  the  bearings  as  well  as  other  parts  of 
the  mechanism. 

ji.  The  ability  to  increase  rpm  allows  the  rotational  system  to  test 
at  high  discharges,  provided  the  fixture  can  handle  the  peak 
forces.  Consequently,  it  is  more  flexible. 

Precise  timing  for  the  start  of  the  test  is  critical  to  getting 
meaningful  results. 

_f.  Neither  system  provides  enough  time  at  high  discharge  to 
thoroughly  test  the  regulator  performance. 
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APPENDIX  A:  ROTATIONAL  TEST  RESULTS 


1.  The  figures  presented  in  this  appendix  show  the  rotational  testing 
system's  performance  at  the  design  shaft  speed  of  424  rpm.  Other  analysis 
(not  presented)  showed  that  the  design  discharge  (26.4  gps)  would  be  attained 
only  when  pumping  action  starts  at  a  shaft  angle  of  ninety  degrees.  This 
analysis  shows  the  system's  performance  under  those  conditions. 

2.  Figure  A1  provides  a  tabulation  of  various  system  parameters  versus 
time.  Figures  A2  through  A6  are  plots  of  those  same  parameters  versus  time. 
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APPENDIX  B:  ROTATIONAL  TEST  RESULTS,  PUMP  START  AT  ZERO  DEGREES 


1.  The  figures  presented  in  this  appendix  show  the  rotational  testing 
system's  performance  at  two  shaft  speeds.  Other  analysis  (not  presented) 
shows  that  the  design  discharge  (26.4  gps)  would  be  attained  at  the  design  rpm 
(424)  only  when  pumping  action  started  at  a  shaft  angle  of  ninety  degrees. 
However,  at  that  point  the  piston  and  fluid  are  moving  at  the  maximum  velocity 
making  pump  starting  very  difficult  under  those  conditions.  This  analysis 
shows  two  sets  of  results  for  starting  the  test  at  zero  degrees. 

2.  Figure  B1  provides  a  tabulation  of  various  system  parameters  versus 
time.  That  output  shows  the  peak  discharge  reduced  to  about  16  gps  when  the 
shaft  speed  at  pump  start  is  424  rpm.  Figures  B2  through  B6  show  plots  of 
discharge  and  other  system  parameters  versus  time  for  the  424  rpm  test. 

3.  Figures  B7  and  B8  show  the  peak  system  discharge  to  be  at  the  design 
discharge  when  the  shaft  speed  at  pump  start  is  increased  to  577  rpm.  Figures 
B9  through  B1 2  show  other  system  parameters  versus  time  at  577  rpm. 
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388.2 

40103. 

12.79 

-26317 

4  00 

0.0178 

42 . 92 

362.2 

34907. 

15.02 

-33101 

500 

0 . 0222 

52.16 

330. B 

29121 . 

15.98 

-38542 

600 

0.0266 

60.50 

295.3 

23199. 

15.78 

-42654 

70  0 

0.0311 

67.86 

256.7 

17529. 

14. 66 

-45569 

800 

0 . 0355 

74.15 

216.0 

12409 . 

12.85 

-47493 

900 

0.0400 

79.34 

173.9 

8047, 

10,61 

-48662 

1000 

0.0444 

83.40 

131.1 

4572. 

8. 10 

-49304 

1100 

0.0438 

86.31 

87.9 

2053. 

5 , 46 

-49614 

1200 

0.0533 

38.07 

44.4 

525. 

2.77 

-4°739 

1300 

0 . 0577 

88 . 66 

1 .0 

0. 

0.06 

-49771 

PUMPING  ENDS  AT  0.058  SEC. 

PUMP  CUTOUT  ANGLE  =  88.7  DEGREES 

NUMBER  OF  TIME  STEPS  -  1303 
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FIGURE  B5 


tine  in  niu-sec 


TEST  START  ANGLE  ?  RPM  »? 

-0  577 


N 

T( SEC) 

THETA 

RPM 

KE  (FT-LB) 

Q  (GAL/SEC) 

TORG  <  FT-LB ) 

0 

0  . 

0  . 

577.0 

88582. 

0.06 

0  . 

100 

0 . 0060 

20.79 

566 . 8 

85472  . 

12.12 

-17136. 

200 

0.0121 

40 . 86 

537 . 4 

76351  . 

21 .36 

-31771  . 

300 

0.0131 

59.57 

493.0 

64679  . 

26 . 09 

-42233. 

400 

0.0242 

76.47 

438.9 

51265 . 

26.44 

-48063. 

500 

0.0302 

91 . 32 

380 . 4 

38501  . 

23.78 

-49849 . 

600 

0.0363 

104.03 

321 .6 

27511  . 

19.65 

-48728. 

700 

0.0423 

114.64 

265.2 

18707. 

15.27 

-45915. 

800 

0.0483 

123.28 

212.6 

12025. 

11.31 

-42414 . 

900 

0.0544 

130.09 

164.2 

7177. 

8.02 

-38934 . 

1000 

0.0604 

135.22 

119.8 

3817  . 

5.40 

-35925. 

1100 

0.0665 

133.80 

78.5 

1639. 

3.31 

-33642. 

1200 

0 . 0725 

140.92 

39.4 

413. 

1 .59 

-32219 . 

1300 

0.0736 

141.66 

1 .5 

1  . 

<1 

o 

o 

-31718. 

PUMPING  ENDS 

AT 

0.079  SEC 

♦ 

PUMP  CUTOUT  ANGLE  =  141.7  DEGREES 

NUMBER  OF  TIME  STEPS  =  1304 
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APPENDIX  C:  ROTATIONAL  SYSTEM  PROGRAM  "JPUMP 


1.  This  appendix  presents  a  listing  of  program  JPUMP  that  was  used  for 
the  idealized  dynamic  analysis  of  the  rotational  test  system.  This  program 
produced  the  tabulations  and  plots  presented  in  Appendixes  A  and  B. 


JF’UMF' 

IOCS  TITLE  FILED  IN  ROSSCOBRA/ JPUMP 

20C  F'UMP  RUN-TIME  TEST  PROGRAM 

30C  JAY  CHEEK ,  BHD,  SL,  WES,  14  JAN  1986 

TOC 

50 C  HP  PISTON  DAMETER  (IN.) 

60C  US  PUMP  SHAFT  DIAMETER ( IN . ) 

7 0 C  AP  EFFECTIVE  AREA  OF  PISTON  (FT**2) 

530C  SP  STROKE  OF  THE  PISTON  (IN) 

90C  DF  FLYWHEEL  DIAMETER  (IN.) 

100C  WF  FLYWHEEL  WEIGHT  (LBS.) 

110C  P I  3.1415927 

1 20C  DTR  FACTOR  TO  CONVERT  DEGREES  TO  RADIANS 

1 30C  G  32.2 

1 4 OC  P  WORKING  PRESSURE  IN  THE  PUMP  (P/IN#*2) 

1 50C  FM  FLYWHEEL  MASS 

1 60C  F.'MI  ROTATIONAL  MOMENT  OF  THE  FLYWHEEL 

1 7 0 C  RPM  INITIAL  REV.  PER  MIN.  OF  THE  FLYWHEEL 

1 8 0 C  AVF  CURRENT  ANGULAR  VELOCITY  OF  THE  FLYWHEEL  (RAD./ 

1 90C  AVF I  INITIAL  ANGULAR  VELOCITY  OF  THE  FLYWHEEL  (RAD,/ 

2 0 0 C  EK  KINETIC  ENERGY  OF  THE  FLYWHEEL 

2 1 OC  THETAI  SHAFT  ANGLE  WHEN  PUMPING  TEST  STARTS  (DEG) 

220C  THETA  CURRENT  SHAFT  ANGLE  DURING  PUMP  TEST 
230C  ZMI  INITIAL  MOMENTUM  OF  THE  FLYWHEEL 

2 4 0 C  ZM  CURRENT  MOMENTUM  OF  THE  FLYWHEL 

2 5 0 C  ZIMP  TOTAL  IMPULSE  ON  FLYWHEEL  FROM  PUMP  ACTION 

2 6 0 C  TDEL  TIME  STEP (SEC) 

270C  T I  CURRENT  TIME  SINCE  START  OF  THE  TEST. 

28 OC  CL  CONNECTING  ROD  LENGTH. 

290C  TOR  TORQUE  APPLIED  DURING  ONE  TIME  STEP 

3 0 0 C  ALPHA  ANGULAR  DECELERATION  DURING  ONE  TIME  STEP, 

3 1 OC  Q  DISCHARGE  OF  THE  PUMP  IN  GAL/SEC. 

320C  XP  D 1 3 T .  FROM  CL  OF  CRANK  TO  PUMP  AT  STEP  N 

330C  XFL  DIST.  FROM  CL  OF  CRANK  TO  PUMP  AT  STEP  N-l 

3  4  0C 
35  OC 

360  DIMENSION  V( 1500,6) 

370  DATA  PI/3.1415927/,  G/32.2/,  DP/16.0/,  DS/9.0/, 

380  2  DTR/57. 295780/,  WF/1000.0/,  SP/2.0/,  P/4350.0/, 

390  l  DF/30.0/,  CL/30.0/ 

4  0  0  C 

4 1 OC  PLOT  SETUP 

420  CALL  USTART 

430  CALL  UPSET ( 'SPEED' ,  120.) 

440  CALL  U3ET ( 'SMALL  '  ) 

450  CALL  IJSET  (  '  XBOTHL ABELS  '  ) 

4  60  CALL  USET ( ' YBOTHL ABELS  '  ) 

470  CALL  UERASE 

480  CALL  UHOME 

490  CALL  UhLF'HA 

5  0  0  C 


C2 


GO  to 


■.nr  im  '.m  im  m  'J*J  "VliV'.ll  l  ■v 1  »v  **T  v  V  r"  .  »  -  *  .  '  '■'  ■ 
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JPUMP 

510C-  GET  INITIAL  CONDITIONS 

02  0  100  WRITE  (or  110)  '  TEST  START  ANGLE  ,  RPM  ' 

530  110  FORMAT  (0) 

340  READ  (5»110)  THETA  I  >  RPM 

3 5 0 C  PISTON  AREA 

360  AP  -  PI  *  (DP  *  DP  -  DS  *  DS)  /  576.0 

570C  MASS  OF  THE  TWO  FLYWHEELS 

580  FM  -UF  *  2.0  /  G 

590C  ROTATIONAL  MOMENT  OF  INERTIA 

600  RMI  -  FM  *  DF  *  DF  /  1152.0 

6  IOC  INITIAL  ANGULAR  VELOCITY 

620  A OF  I  =  RPM  *  PI  /  30.0 

63  OC 

6 4 0 C  CALC  CONSTANTS  FOR  THE  INTEGRATION  LOOP. 

6 50C  CRANK  ARM  LENGTH  IN  FT 
660  TH  =  SP  /  24.0 

6 7 OC  SQUARE  OF  THE  CONNECTING  ROD  LENGTH 
630  CLS  =  CL  *  CL  /  144.0 

690C  TIME  STEP  IN  SEC. 

700  TDEL  =  AVF  I  /  1000000.0 

7  IOC  PUMP  FORCE  TIMES  CRANK  ARM  LENGTH 

720  FR  =  AP  *  P  #  TH  *  144.0 

730C  STEP  COUNTER 
740  N  =  0 

7 5 0 C  INITIAL  ANGULAR  VELOCITY  OF  THE  CRANKSHAFT  (RADIANS/SEC) 


7  6  0 

A  V  F  -  A  V  F  I. 

7  7  0  C 

CRANKSHAFT  ANGLE  AT  START  OF  PUMPING 

(0  DEG 

IS 

7  8  0  C 

TOP  DEAD  CENTER  OF  STROKE) 

79  0 

THETA  -  THETAI 

3  0  0  C 

NUMBER  OF  POINTS  TO  BE  PLOTTED 

8  1  0 

NVP  =  0 

S20C 

330C 

HEADING . 

‘ 

840 

WRITE  (6*120) 

Q  17" 

U  J  J 

120  FORMAT  ('  N  T(SEC)  THETA 

RPM 

KE  (FT-LB) 

Q  ' 

- 

360 

S  *  '(GAL/SEC)  TORQ  (FT-LB)’) 

8  7  0  C 

8S0C 

START  THE  PUMP  TEST  LOOP 

890C 

ELAPSED  TIME  IN  SEC 

900 

130  T I  =  N  *  TDEL 

9  IOC 

CURRENT  CRANKSHAFT  ANGLE  IN  RADIANS 

920 

TR  =  THETA  /  DTR 

930C 

KINETIC  ENERGY  IN  FT  LBS 

940 

EK  =  RMI  *  A VF  *  A VF  /  2.0 

9  50C 

TEST  FOR  SPECIAL  CALCULATION  FOR  Q  AT 

START 

960 

IF  ( N  .HE.  0)  GO  TO  140 

9  7  0  C 

CALC  FUMP  POSITION  JUST  BEFORE  THE  FIRST  STEP. 

980 

TRR  =  TR  --  A  V  F  *  TDEL 

9  90 

ST  -  SIN (TRR) 

1000 

TEMPI  -  TH  *  COS (TRR) 

C3 


JPUHP 

1010  TEMP2  =  SGRT'CLS-<TH*ST)**2) 

1020  XPL  --  TEMP 2  --  TEMPI 

1030  140  ST  =  SIM(TR) 

1040C  X  COMPONENT  OF  THE  CRANK  ARM 
1050  TEMPI  =  TH  *  COS(TR) 

1060C  PROJECTION  OF  CONN*  ROD  ON  X  AXIS 
.1.070  TEMP2  =  SQRT(CLS-(TH*ST>**2> 

1 0 8 0 C  POSITION  OF  THE  PUMP  AT  STEP  N 
1090  XP  TEMP 2  -  TEMPI 

1100C  PUMP  DISCHARGE  DURING  THIS  STEP 
1110  Q  =  ABS(XP-XPL)  *  AP  *  7.5  /  TDEL 

1120  XPL  =  XP 

113QC  SLOWDOWN  TORQUE  ON  CRANK  DUE  TO  PUMPING  FORCE 
1140  TOR  =  --  ABS<FR#ST*(1 . 0-TEMP1/TEMP2 )  ) 

1150C  ANGULAR  DECELERATION 
1160  ALPHA  =  TOR  /  RMI 

1170  IF  ( MOD ( N , 1 00 )  ,EQ.  0)  WRITE  (6,150)  N  ,  TI,  THETA,  AVF  / 

1130  &  2.  /  PI  *  60,  EK ,  Q ,  TOR 

1190  150  FORMAT  (IX?  14.*  F  3 . 4  ,  F3.2,  FB.l,  FI  0.0,  F 1 4  *  2  *  FI  4.0) 

1200C 


1210  C 

SAVE  EVERY  FIFTH  POINT  FOR 

PLOTTING 

12  2  0 

IF  ( M  0  D ( N , 

5)  .NE.  0)  GO 

TO  160 

12  30 

NVP  =  NVP 

+  1 

1  2  4  0 

V  (  N  V  P  ,  1  )  - 

ri  *  iooo.o 

1 2  5  0 

V ( NVP , 2 )  = 

a 

1260 

V ( NVP , 3 )  - 

TOR 

1 2  7  0 

V  (  N  V  P  ,  4  )  -- 

EK 

1230 

V  (  N  V  P  ,  5  )  ~ 

THETA 

.1.29  0 

V < NVP, 6  )  - 

AVF  /  2.0  / 

PI  *  60.0 

130 0 C  NEW  ANGULAR  VELOCITY  DUE  TO  SLOWDOWN  TORQUE 

1310  160  A V F  -  A V F  +  ALPHA  *  TDEL 

1320C  COUNT  STEF 

1330  N  -  N  +  1 

1340C  STILL  ROTATING? 

1350  IF  ( A  V F  . LT.  0.0)  GO  TO  170 

1360C  YES.  ADVANCE  TO  NEXT  CRANK  ANGLE 
1370  THETA  -  THETA  +  A VF  *  TDEL  *  DTR 

1330  GO  TO  130 

1390C 

1 4  0  0  C  DONE 

1410  170  WRITE  (6,180)  TI,  THETA  -  THETA  I ,  N 

1420  130  FORMAT  ('  PUMPING  ENDS  AT  ' ,  F10.3,  '  SEC.'  / 

.1430  S  '  PUMP  CUTOUT  ANGLE  =  '»  F8.1,  '  DEGREES'  / 

1440  S  '  NUMBER  OF  TIME  STEPS  =  .15  ////  ) 

1450  CALL  UFAUSE 

1460C 

1470C  PLOT  THE  CURVES 

1430  CALL  PLOTPU ( V »  NVP ,  DP,  DS,  P,  DF-  WF,  SP«  CL,  RPM,  THETA  I  ) 

1490  GO  TO  100 

1500  END 


C4 


JPUMP 

1510  SUBROUTINE  PLOTF'UCV?  NV?  DP?  DS?  P-  DF  ?  WF?  SP ?  CL?  R  P  M  ? 

1520  l  FHETAI) 

.1 5  3  0  C 

1540C  PLOT  ANY  OF  THE  5  ARRAYS  AS  A  FUNCTION  OF  TIME  (ARRAYl), 

1 5 5 0 C  JAY  CHEEK  ?  S H D ?  3 L  r  WES  »  IS  JAN  1 9 8 S 
1  5  6  0  C 

1570  CHARACTER*20  TC20) 

1580  DIMENSION  VC 1500,6) 

1590  DATA  T  C  1  )  /  '  T I M  E  IN  MILI-SEC\  '/?  T(2) 

1600  S  /"Q  IN  GAL  PER  SEC\  '/?  TC3) 

1610  l  /'TORQUE  IN  FT  LB\  '/?  T(4> 

1620  S  /'KE  IN  FT  LB\  '/?  T(5> 

1630  l  /"CRANK  ANGLE  IN  DEG \  '/?  TC6) 

1640  l  ./'REV  PER  MINX  '/ 

1650C  START  PLOTTING 
1660  FN  =  NV 

1670  CALL  UPSET ( ' XLABEL  '  ?  T(l>) 

1680  DO  120  I  =  2?  6 

169  0  C 

17  00 C  IS  THIS  GRAPH  NEEDED? 

1710  WRITE  (6? 100)  TCI) 

1720  100  FORMAT  C'  PLOT  '?  A20?  '?') 

1730  CALL  IAN SRC  1ST) 

1740  IF  (1ST  .EQ.  0)  GO  TO  120 

173  0  C 

1760C  YES.  DO  IT, 

1770  CALL  UPSET C 'YLABEL ' ?  TCI)) 

1730  CALL  UERASE 

1790  CALL  UDAREAC4.0?  14,0.-  0,0.  10.0) 

1800  CALL  U3ETC 'GRIDAXIS' > 

1810  CALL  UPLQTKVC1?  1),  MCI,  I),  FN) 

1820  CALL  UHOME 

1330  CALL  UALPHA 

1840  WRITE  (  6  ?  1 1 0  )  DP?  DS  ?  P  ?  DF  t  WF  *  2,0.  SP  /  2.0?  CL?  RPM? 

1350  l  THETAI 

I860  110  FORMAT  C'  PROGRAM  ROSSCOBRA/ JPUMP  --  JAY  CHEEK?  SMD  ?  ' 

1870  l  7  'SL?  WES?  16  JAN  1986'  /  '  PISTON  DIAMETER  =  '? 

1830  l  F3.2?  '  IN'  /  '  SHAFT  DIAMETER  ='?  F8.2?  '  IN'  / 

1390  S  '  PUMP  PRESSURE  ='?  F8.2?  '  PS  I'  / 

1900  l  '  FLYWHEEL  DIAMETER  = ' ?  F 8 . 2 ?  '  I N '  / 

1910  S  '  FLYWHEEL  WEIGHT  ='?  F3.2?  '  LB'  / 

1920  3  '  CRANK  THROW  = ' ?  F6.2?  '  IN'  /  '  CONNECTING  ROD  LENGTI 

1930  l  ?  F3.2?  '  IN'  /  '  INITIAL  RPM 

1940  i  ?  F  8 , 2  /  '  CRANK  ANGLE  AT  START  ='?  F8.2?  '  DEG') 

1950  CALL  UPAUSE 

I960  120  CONTINUE 

1970  RETURN 

1980  END 


0-  LI 


'  U  M  P 

.1  990  SUBROUT  I  NE  IANS R  <  I  WHAT  ) 

2 0 0 0 C  SEE  WHETHER  THE  USER  GIVES  A  Y  (OR  A  CARRIAGE  RETURN) 

20 IOC  FOR  YES  OR  A  N  FOR  NO  TO  A  PREVIOUSLY  ASKED  QUESTION. 

'7  C i  '7  Cf  0 

2  0  3 0  C  CODE  FILED  IN  R0S3C0BRA/ JHEST-S . 

2 0 4 0 C  JAY  CHEEK  t  SMD,  3L ,  WES  5  DEC  1981 
2050  C 

2060  DATA  I E  L  K /  '  '/»  IYES/'Y  '/,  NO/'N  '/ 

2070  .100  I  WHAT  =  0 

2080  READ  (5»110)  II 

2090  110  FORMAT  (A4) 

2100  IF  (II  .EG.  NO)  RETURN 

2110  I  WHAT  =  1 

2120  IF  (II  ,  EQ.  I  YES)  RETURN 

2130  IF  (II  .EQ.  IBLK)  RETURN 

2140  WRITE  ( 6  > 1 2 0 ) 

0  120  FORMAT  ('  ERROR!  ONLY  Y  OR  RETURN  (FOR  YES)  OR  N  (FOR 

0  £  r  'NO)  ALLOWED*  RETRY') 

2170  GO  10  100 

2130  END 


C6 


APPENDIX  D:  METHOD  FOR  CALCULATING  SLOWDOWN  OF  FLYWHEEL 
Initial  Conditions 

1.  As  shown  in  Figure  D1 ,  a  flywheel  (solid  disc)  of  weight  W  is 
rotating  at  an  initial  rpm.  The  flywheel  diameter  is  D.  The  flywheel  is 
connected  to  a  crankshaft  whose  offset  (throw)  is  T.  A  connecting  rod  of 
length  C  connects  the  crank  to  the  pump  that  resists  motion  with  a  constant 
force  F.  That  force  is  directed  on  a  line  from  the  crankshaft  center-line 
through  the  center-line  of  the  pump  shaft. 


Pump's  End  Plates  and  Valves  not  shown 


Idealized  Pump  and  Drive  System 
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Calculating  Proc'jdure  For  Flywheel  Slowdown 

2.  The  following  steps  are  used  to  calculate  the  flywheel  slowdown: 
a.  Calculate  the  Rotational  Moment  of  Inertia  (I). 


i  -  *  (?)  1 

g  2  2 


b.  Calculate  the  Angular  Velocity  (w). 
(rpm) 


a)  =  2it 


60 


and  choose  a  small  time  step  (At)  so  that  the  wheel  will  rotate 
less  than  .1°  at  the  initial  rpm  during  time  At. 


£.  At  each  time  Step  i,  Calculate: 

1.  The  angle  traversed  (A8)during  this  time  step  (constant 
angular  velocity  (u>)  is  assumed). 


A0  =  w.At 
l 


2.  The  current  angular  location  (©)  of  the  crank. 


0.  =  0._+  A0 


3.  The  component  of  F,  acting  through  the  connecting  rod 


i 


i 


6.  The  new  angular  velocity  due  to  the  slowdown  torque  from 
the  pumping  force. 


a).  =oj.  f  a  At 
l  1-1 


7.  Add  1  to  i  and  repeat  steps  Cl  through  C7  until  uk  £  0. 

8.  Done. 

Calculating  Moment  Applied  To  The  Crankshaft 

3.  As  illustrated  in  Figure  D2:  T  is  crankshaft  throw,  F  is  force  to 
operate  the  pump,  C  is  connecting  rod  length,  E  is  the  moment  arm, 

0  is  crankshaft  angle,  and  G  is  the  force  normal  to  E. 


Calculating  the  Torque  (r)  about  point  Q. 


T  =  G  E 
E  =  A  sin 
A  =  B  -  Tv 


B  =  (C*  -  Tyd) 
Ty  =  T  sin  0 
Tx  =  T  cos  0 


G  = 


COS  p 

B 


r  ■  t  <B  -  V  c  ■  F  V  --r> 


COS  p  =  - 


- 0' 


sin  p  =  -£*- 


T  sin© 


J 

* 

5 

■j 


r  =  F  T  sin  0  [  1  -  - 


T  cos  0 


(C2-  T2  sin20)2 


]  for  C  >  T 


D3 


L'- V- 


*v 

•V*V-*aV*V*  -  • . 
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APPENDIX  E:  ANALYSIS  OF  THE  DROP  TEST  SYSTEM 


1.  This  appendix  contains  the  output  of  the  drop  test  analysis  program, 
JPEND  (Figure  El),  the  listing  of  the  program  (E3),  and  the  development  of  the 
equations  used  (E4  and  E5).  The  output  is  given  for  both  1,000  and  2,000 
pound  balls  impacting  on  the  pistons  of  various  diameter  cylinders.  The  drop 
height  is  paired  with  the  cylinder  diameter  to  produce  the  rated  discharge 
peak  of  26.4  gps.  Discharge  pressure  is  assumed  to  be  constant  at  4,350 
psi.  Note  that  an  increase  in  ball  weight,  holding  other  parameters  constant, 
serves  to  increase  the  test  time  (t)  directly  as  the  ratio  of  the  new  weight 
to  the  old 

WN 

tM  =  — —  t  where  subscripts  0  and  N 

Wq  O 

refer  to  old  and  new  values. 

Figure  E2  gives  the  discharge  versus  time  curves  for  several  ball  weight, 
cylinder  diameter  and  drop  height  combinations. 


Output  of  the  Drop  Test  Ana lysis  Program  ' 'JPEND 


DROP  TEST  OF  0 ALOE 

FOR  BALL  LIE  I GHT  OF  1000.  LB  PEAK  DISCHARGE  OF  26.4  GAL  /  SEC 
AND  DISCHARGE  PRESSURE  OF  4350.  PSI 


CYLINDER 

DROP 

RETARD 

IMPACT 

DISCHARGE 

KINETIC 

DIAMETER 

HEIGHT 

FORCE 

VELOCITY 

TIME 

ENERGY 

(INCHES) 

(  FEET  ) 

(POUNDS) 

(FT  /  SEC) 

(SECONDS) 

(FT  -  LBS) 

3.50 

4  3. 10 

41352 . 

52 . 68 

0 . 039 

43101 . 

3,75 

32 .71 

48044 . 

45.89 

0.030 

32706 . 

4.00 

25.26 

54664 . 

40.34 

0.023 

25265. 

4.25 

19.82 

61710. 

35.73 

0.013 

19S24  . 

4,50 

15,77 

69134 . 

31 .37 

0.014 

15773. 

4,75 

12.71 

7 7084 » 

28,60 

0.012 

12705. 

5.00 

10.35 

35412, 

25.82 

0.009 

10348. 

DROP  TEST  OF  VALUE 

FOR  BALL  WEIGHT  OF  2000.  LB  PEAK  DTSCHARGE  OF  26.4  GAL  /  SEC 
AND  DISCHARGE  PRESSURE  OF  4350.  PSI 


CYLINDER 

DROP 

RETARD 

IMPACT 

DISCHARGE 

KINETIC 

DIAMETER 

HEIGHT 

FORCE 

VELOCITY 

TIME 

ENERGY 

(INCHES) 

(FEET  ) 

(POUNDS) 

(FT  /  SEC) 

(SECONDS) 

(FT  -  LBS) 

3 . 50 

4  3.10 

41852 . 

52,68 

0,078 

86201 . 

3 , 75 

32.71 

48044 , 

45.89 

0,059 

65413 . 

4 , 00 

25 . 26 

54664 . 

40,34 

0.046 

50530 . 

4.25 

19.82 

61710. 

35.73 

0 . 036 

39649  . 

4 , 5  0 

15 . 77 

69184 . 

31,87 

0,029 

31545  . 

4 . 75 

12.71 

77084  . 

28.60 

0,023 

25410  . 

5,00 

10.35 

85412  . 

25 , 82 

0.019 

20697. 

FIGURE  E2 


PEND 

10C$TITLE  FILED  IN  ROSSCOBRA/ JPEND 

20C  CALC  OF  PENDULUM  TEST  OF  VALUE  PERFORMANCE, 

30 C  JAY  CHEEK  SMD,  SL ,  UES,  17  JAN  1986 
40C 

50  DATA  PI/3. 1415926/,  P/4350. 0/»  0/26.4/ 

60  C 

70  DO  130  J  =  1»  2 

SOC  WEIGHT  OF  BALL 
90  WT  ----  J  *  1000.0 

100C  DIAMETER  OF  HYDRAULIC  CYLINDER  (ONE  STROKE  PUMP) 

110  D  -  3.250 

120  WRITE  <  6 » 100)  WT,  Of  P 

130  100  FORMAT  (  ////  '  DROP  TEST  OF  VALVE'  / 

140  8  '  FOR  BALL  WEIGHT  OF ' >  F  6 . 0  , 

150  8  '  LB  PEAK  DISCHARGE  OF',  F5.1,  '  GAL  /  SEC'  / 

160  8  '  AND  DISCHARGE  PRESSURE  OF',  F6.0,  '  PSI '  // 

170  8  '  CYLINDER  DROP  RETARD  IMPACT 

180  8  ,  'DISCHARGE  KINETIC'  / 

190  8  'DIAMETER  HEIGHT  FORCE  VELOCITY  TIME  ' 

200  8  ,  '  ENERGY'  / 

210  8  '  (INCHES)  (FEET)  (POUNDS)  (FT  /  SEC)  ' 

220  8  ,  '(SECONDS)  (FT  -  LBS)') 

23  OC 

2 4 0 C  CALC  FOR  SEVERAL  CYL .  DIAMETERS 
250  DO  120  1=1,7 

2  6  0  C  DIAMETER 

270  D I  -  D  +  I  *  .25 

28 OC  PISTON  (CYLINDER)  AREA 

290  A  =  PI  *  D I  *  D I  /  4,0 

3 0 0 C  MASS  OF  THE  BALL 

310  ZMASS  =  WT  /  32.2 

3 2 0 C  DROP  HEIGHT  TO  PRODUCE  PEAK  DISCHARGE 

330  H  =  9.230  *  Q  *  Q  /  DI  **  4 

34 OC  IMPACT  VELOCITY 

350  V  =  SORT (2. 0*32. 2*H) 

360C  RETARDING  FORCE  OF  THE  CYLINDER 
370  F  =  P  *  A 

33 OC  TIME  OF  DISCHARGE 
390  T I  =  ZMASS  *  V  /  F 

400C  KINETIC  ENERGY 

410  EK  =  ZMASS  *  V*  V  /  2.0 

4  2  0  C 

4  3  0  C  RESULTS 

440  WRITE  (6,110)  DI ,  H,  F,  V,  TI ,  EK 

450  110  FORMAT  (IX,  F7.2,  F10.2,  F12.0,  F10.2,  F10.3,  F10.0) 

460  120  CONTINUE 

470  130  CONTINUE 

480  STOP 

490  END 


E3 


Equation  For  Drop  Height 

2.  Calculating  combinations  of  cylinder  diameter  (D)  and  drop  height 
(H)  that  yield  a  specific  value  of  fluid  discharge  Q. 

1 

a_.  Impact  velocity  =  Uj  =  (2gH)  2 


where  U  is  in  ft/sec 
H  is  in  ft 
g  is  in  ft/sec* 


V 

_b.  Discharge  =  Q  =  — —  where  V  is  volume  in  ft. 


Q  -  1 -  AU. 

At  At  I 


where  A  is  cylinder  area  =  jj  D2 
D  is  cylinder  diameter  in  ft 


2g  (jj  D2)2 


1 .23^gD 


where  Qj  is  discharge  at  impact,  i.e.,  peak  discharge, 


d.  Converting  for  Q  in  gps,  and  D  in  inches. 


4> 

D  4 

1.234  (yj)  g 


9.280  Q. 


Equation  For  Total  Discharge  Time 

3.  As  shown  on  E4,  Drop  Weight  (H)  relates  to  Peak  Discharge  (Qj)  and 
cylinder  diameter  (D)  by: 


„  9.280Q, 

“7"1 


Where  Q  is  in  gps,  D  is  in  inches, 


Conservation  of  momentum  gives  MUj  =  Ft 


Where  t  is  the  time  in  seconds  that  the 
constant  force  F  exists,  M  is  the  ball 

mass,  W  is  the  ball  weight  and  M  =  W. 

g 


F  is  the  constant  piston  force  = 


2 

IT  D  P 


Where  P  is  pressure  in  psi,  F  is  in  lbs. 


U j  =  impact  velocity  =  (2gH)  2 


Total  discharge  time  (t)  is 


W  ; 

g  (2gH)‘ 


2g9.280Q" 


!d2P 


-  .9667 


Where:  W  is  weight  in  lb,  D  is  diameter  in  inches,  Q  is  discharge  in 
gps,  P  is  pressure  in  psi,  and  t  is  time  in  sec. 

For  P  =  4,350  psi,  Qj  =  26.4  gps 


t  =  .005866  W 


E5 


;■!]  ryyyy. 
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